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Abstract 

In this paper we introduce the data from mineral water probe with errors in both variables. 
For this case we apply our orthonormal polynomial expansion(OPEM) method to describe the 
data in the new error corridor. It received the approximating curves and their derivatives 
including the errors by weighting approach. The numerical method and approximation results 
are presented and discussed. The special criteria are carried out for orthonormal and evaluated 
from it usual expansions. The numerical results are shown in tables and figures. 
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1 Introduction 



The water spectra method applies a drop taken from a water probe to measure the water's 
state spectrum. In the special experiment the drop is placed on a hostaphan folio- Figure 1[1]. 
JS^ ! During the whole process of evaporation of the drop, one measures at equal time intervals the 

^ drop contact angle with the folio. On X-axis one has the values of the contact angles within 

fixed angular intervals and on Y-axis the frequency of measurements of these angles. 

To compare different state spectra one normalizes each spectrum dividing its Y-values by 
the number of all measurements and thus obtains a probability distribution . One can change 
the function (f){6) on the independent angle 6 to the function of energy variable F(E) using the 
following Antonov transformation [2]: 

f{E) = 60(^)/v'l-(l + 6E)2, 

where 

b = i{i + cosieo))h- 

Here I = 5.03.1018m~^ is the density of water molecules in the surface layer, 7 is the surface 
tension, ^0 - the initial contact angle. 

The so obtained graph after measurements by method in [1] is referred to as energy spectrum 
F{E) of the probe. E is the energy of Hydrogen bond of investigated water. 

On Figure 2 the dependent variable contains the values of the water Hydrogen bond energy. 
Here we present new detailed information about given data and their errors on both variables 
of water probe. 

The method of water spectra is sensitive to treatment by physical fields as 7-ray treatment 
of water [1, 2] and to environmental changes of the ecosystem on different water probe [3, 4, 5, 
6, 7, 8]. In the present paper we approximate another natural water data taken from a water 
spring in Bulgaria near the village Lenovo. 
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Figure 1: Experimental setup 




Figure 2: Experimental data from (B) with their errors 

2 Main problem definition 

• To find the best approximation curves of measured water data on Fig. 2 taking into account 



• To extend our original Orthonormal polynomial expansion method (OPEM), according 
some criteria, to evaluate orthonormal description of given data; 



• To find the best approximating curve with usual polynomials, evaluated by orthonormal, 
according some criteria. 

3 Numerical method— OPEM "total variance" 

Let the {Ei,Fi,i — 1, ...M} are arbitrary pairs of monitoring data E = Ei and F = Fi, 
introduced in section 2. They are given with experimental errors in both variables-o'(Fj) and 
a{Ei). Consider the total uncertainty (total variance) S^{E,F) [9, 10, 11], associated with 
(E, F) 

SI = a\F,) + (1) 

according the ideas of Bevington (1977) [9], where his proposal is to combine the errors in both 
variables and assign them to dependent variable. One defines the errors corridor C{E,F), 
which is the set of all intervals 

[F(E)-S(E,F),F{E) + S(E,F)], (2) 



3.1 orthonormal expansion criteria 

The first criterion to be satisfied, is that the fitting curve should pass within the errors corridor 
C{E,F). In the cases of errors only in F, (i.e. cr(£') = 0,a{F) ^ (0)) the errors corridor 
C{E,F) reduces to the known set of intervals 

[F-<j{F),F + a{F)l (3) 

for any F. The second criterion is, that the fitting curve F"''^{Ei) satisfies the expression 

M 

= - F{E,)fl{M - L),w, = 1/Sl (4) 

i=l 

should be minimal (L-number of polynomials). The preference is given to the first criterion. 
When it is satisfied, the search of the minimal chi- squared stops. Some details of the calcula- 
tion procedure are given in Forsythe's paper [12] and in our works [13, 14, 15]. 

Our procedure gives results for approximating function by two expansions : of orthogonal 
coefficients {aj} and usual ones {cj} with optimal degree L: 



,(m) 

i 

i=0 i=0 



The orthogonal coefficients are evaluated by the given values Fj, weights and orthogonal poly- 
nomials: 

M 

a,^Y.FkWuPt\Ek). (6) 



Our recurrence relation for generating orthonormal polynomials and their derivatives (m = 
1, 2...)( or their integrals with m=-l,-2,-3,...) are carried out by: 

P^TiiE) = 7^+l[E - ^^+.)Pt\E) - (1 - 5,o)v,PthE) + mPt-'^E)i (7) 

where jii and z/j arc recurrence coefficients, and 7j is a normalizing coefficient, defined by scalar 
products of given data. One can generate P"^{E) recursively. The polynomials satisfy the 
following orthogonality relations: J2fiiWiPj.^\Ei)Pi'^\Ei) = 6k,i over the discrete point set 
{Ei, i = 1, 2, . . .}. All the calculations for the sake of uniformity are carried out for E in[-l,l], 
i.e. after the input interval is transformed to the unit interval. We remark some advantages 
of OPEM: It uses unchanged the coefficients of the lower-order polynomials; it avoids the 
procedure of inversion of the coefficient matrix to obtain the solution, the minimal chi-squared 
stops. All these features shorten the computing time and assure the optimal solution by the 
criteria(2) and (4). 

3.2 the usual expansion criteria 

The inherited errors in usual coefficients are given by the inherited errors in orthogonal coeffi- 
cients: 

AQ = (E(cS'=))^)^/^Aa., (8) 

k=l 

M 

Aa, = PhEk)MFk - FriV- (9) 

k=l 

(k) 

where coefficients q are defined by orthonormal expansion of polynomials 

k 

Pk = Y.cf^E\k = Q,..,L (10) 

i=0 

and explicitly constructed by recurrence relation in [13]. 

The procedure is iterative because of the evaluation of derivatives on every iteration step 
and the result of the A;**-th consequent iteration is called below the A;'*-the approximation. We 
note that in every iteration step the algorithm find the best approximation using given before 
criteria. 

We can add the other criteria for optimal number of polynomials for usual expansion. Having 
the La we continue with finding the optimal Lc the minimal value in 

max{ci{L),i = 1, L) (11) 

in usual coefficients through all steps of iterations k^* — 1,2, ..9 or we are asking the mini- 
mal value of the maximal distance between functions, evaluated by orthonormal and usual 
expansions 

max\iF:f^^F:j!")\,k=l,M (12) 
through all iterations. We investigate both criteria, but we prefer the last one. 

4 Approximation results 



The main important results from approximation between 2 10 degrees for iterations 1 -j- 9 are 



Table 1: OPEM approximations results for every step approximation 
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Figure 3: OPEM approximation by 6-th degree orthonormal polynomials (C)(2-nd 
iteration) of experimental water data(B) 

max\{Fa — Fc)\. We see from the Table 1, that from iteration number 2 -=-5 with optimal number 
La = 6 the results are good for both expansions and for usual expansion the 8-th iteration with 
optimal number Lc = 8 they are also good. 

Note: It is very interesting to present on figures the three curves - given(B), approximated by 
orthogonal polynomials(C) and received from it by usual polynomials(D) at different iteration 
steps. 

Below the figures 3,4,5,6 present the different approximations results with 2-nd, 3-rd and 
4-th iterations. 

The Table 2 presents the given and approximating values by OPEM with usual and or- 
thonormal coefficients by calculated optimal degree 5 — th in 8 — th iteration of M = 18 given 
values of following characteristics: energy E , distribution F , aE and ap, and from 5 — th 
column - the approximating values with orthonormal coefficients F^pp^'^^ approximating values 
with usual coefficients F^pp'''\ differences A{Fa,Fc) = {F^^pp'''^ - F^pp"'^), total variance 5(5) 
(equation (1). The Table 2 shows good coincidence between two descriptions. For comparison 
we can see the previous results for OPEM applications in [13, 14, 15, 16, 17]. 




Figure 4: OPEM approximation by 6-th degree(iteration-2-nd) orthonormal poly- 
nomials (C) and received usual expansion(D) of experimental water data(B) 




Figure 5: OPEM approximation by 6-th degree(the 3-rd iteration)-the orthonormal 
polynomials (C) and received usual expansion(D) of experimental water data(B) 




Figure 6: OPEM approximation by 6-th degree(the 4-th iteration) orthonormal 
polynomials (C) and received usual expansion(D) of experimental water data(B) 
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Figure 7: OPEM approximation by 5-th degree(the 8-th iteration) orthonormal 
polynomials (C) and received usual expansion(D) of experimental water data(B) 



Table 



2: OPEM approximation of contact water energy data 
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5 Conclusions 

• We have developed new version of OPEM algorithm and Fortran 77 package to include 
errors in both variables according (2) and (4), defined new '"total variance'" and taking 
into account the respective inherited errors (8) and( 9) in coefficients. 

• The approximating curves are chosen at 2 — nd, 3 — rd, 4 — th approximation step by 
optimal degree La — 6 and at 8-th iteration step by optimal degree Lc — 5 to satisfy the 
proposed criteria (2), (4) and (11), (12). The results show that the orthonormal and usual 
expansions values are close to given ones in the whole interval. 

• Our approximating results with optimal degrees of orthonormal polynomials for contact 
(wetting) angle found by orthogonal and usual coefficients show good accuracy and 
stability, demonstrated from Figures and Tables 1,2. We received suitable descriptions 
of the energy variations useful for further investigations. 

• The presented extended algorithm and package OPEM "total variance" with its accuracy, 
stability and speed can be used in other cases of data analysis (as it it shown in our 
previous papers with earlier versions - for calibration problems in high energy physics 
[18]). 
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